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ABSTRACT: Quantum phenomena such as vacuum polarisation in 
curved spacetime induce interactions between photons and gravity with 
quite striking consequences, including the violation of the strong equiv- 
alence principle and the apparent prediction of 'super luminal' photon 
propagation. These quantum interactions can be encoded in an effec- 
tive action. In this paper, we extend previous results on the effective 
action for QED in curved spacetime due to Barvinsky, Vilkovisky and 
others and present a new, local effective action valid to all orders in 
a derivative expansion, as required for a full analysis of the quantum 
theory of high-frequency photon propagation in gravitational fields. 
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1. Introduction 



In classical electrodynamics in curved spacetime, the interaction of the electromagnetic 
and gravitational fields is encoded in the Maxwell action 



r (0) = -\j dxy/gF^F^ (1.1) 



where the field strength F^ v = D^A V — D^A^ is defined with covariant derivatives. The 
interaction involves only the connection and is independent of the spacetime curvature. 
This action therefore embodies the strong equivalence principle (SEP), which states that 
the laws of physics should be the same in the local inertial frames at each point in spacetime 
and reduce to their special relativistic form at the origin of each of these frames. This is 
assured since the connections (though not of course the curvatures) vanish locally in an 
appropriate frame in a Riemannian manifold. 

The picture is rather different in quantum theory. In quantum electrodynamics in 
curved spacetime, vacuum polarisation loops involving the electron induce interactions 
between photons and the background gravitational field that depend explicitly on the 
curvature, with the necessary scale being set by the Compton wavelength of the electron, 
A r = — . These interactions therefore effectively violate the SEP at the quantum level. 

The simplest effective action which encodes some of the essential physics of vacuum 
polarisation was constructed explicitly by Drummond and Hathrell flj] , who used standard 
Schwinger-de Witt proper time/heat kernel methods (see e.g. as well as an 

independent diagrammatic technique to find 



DH 



r (0) + ^2 / dx V9 (aRF^F^ + bR^F^ x F\ + cR^ Xp F^F Xp 
+ dD^ x D v F v x 



;i.2) 



Here, a, 6, c, d are perturbative coefficients of 0(a), viz. 



1 a 13 a la , la , . 

a = o= c = d = (1-3) 

144 tt 360 tt 360 n 30 n K ' 

As well as being a one-loop result, there are two other significant limitations on Tdh- First, 
it is at most of third order in the generalised 'curvatures' (including both the Riemann 
curvature and the electromagnetic field strengths). Second, it is just the lowest-order term 
in a derivative expansion involving higher orders in 0(— ) acting on the curvatures. If we 
restrict ourselves to the interesting case of terms of 0(F 2 ), which determine the photon 
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propagator and thus the quantum effects on the propagation of light in gravitational fields, 

R a 2 

the first restriction is equivalent to working to first order in O(^), or 0(j%) where L is 
a typical curvature scale, i.e. weak background gravitational fields. In terms of photon 
propagation, the second condition implies a restriction to low frequencies, i.e. lowest order 
in O(-y-) where A is the photon wavelength. 

By far the most dramatic consequence of eq. ( |1 . 2| ) is its prediction of superluminal 
photon velocities. As explained in ref.[|T| (see refs. [[7|,|8|,P, |IT||TT] , |T2| for a selection of subse- 
quent work), the new SEP-violating interactions induce curvature-dependent shifts in the 
light cones, which become 

1 8r 

k 2 A2b + Ac)R^k x + —C^ Xp Wk x a v a p = (1.4) 

where /c M is the photon momentum and is its polarisation. Note that we have written 
eq.( |1.4|) in a convenient form involving the Weyl tensor C^xp. For suitably chosen trajec- 
tories in certain curved spacetimes, this condition remarkably permits k 2 to be spacelike. 

Although photon propagation is the most striking application of the effective action, it 
is clearly important much more generally in encoding the often unintuitive effects of gravity 
on quantum fields. In this paper, we present an extension of the Drummond-Hathrell action 
in which we relax the restriction to lowest order in the derivative expansion, while retaining 
the weak gravitational field approximation. The application of this improved effective 
action to the photon propagation problem and the nature of dispersion in gravitational 
fields is considered separately in ref. |12[ ; here we discuss the technical issues arising in the 
construction of the effective action itself. 

At one-loop order, the QED effective action is given (in Euclidean space) by 

r = r (0) +lndetS(x,a;') (1.5) 

where T(o) is the free Maxwell action and S(x, x') is the Green function of the Dirac 
operator in the background gravitational field, i.e. 

Up -m)S(x,x') = — 6(x,x') (1.6) 

In fact it is more convenient to work with the differential operator corresponding to 
the scalar Green function G(x, x') defined by 

S(x, x') = (ip + m)G(x, x') (1.7) 

so that 

(d 2 + iea^F„ u - ^-R + m 2 )G(x,x') = — —6(x,x') (1.8) 

V 4 / y/g 
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Then we evaluate V from the heat kernel, or proper time, representation 

1 f°° ds 



(o) 



e ~ m s Trg(x,x';s) 



where 



d 

VQ(x, x ; s) = ttG{x, % ; s) 
as 



;i.9) 



(1.10) 



with £/(a;, x'\ 0) = -^=<5(x, a:'). Here, V is the differential operator in eq.(|1.8|) at m = 0. 

Our principal result is summarised in the following expression for the one-loop QED 
effective action, to first order in O 



dx^/g 



4 m z V 

+ Gl RF pv F pv + Gl R p „F» x F v x + G"s R p „ Xp F^F Xp 

+ \(gIrd p f^ x d u f v x 

+ Gl R pv D x F Xp D p F pv + Gl R pv D»F Xp D v F Xp 
+ G* 7 R pv D»D v F Xp F Xp + Gl R pv D»D x F Xp F pv 

+ Gl R pvXp D a F° p D x F pv ) 

(1.11) 



The G n (n > 1) are form factor functions of three operators: 

G n = G 



D 2 



D 2 



D 2 



1.12) 



where the first entry {D 2 ^) acts on the first following term (the curvature), etc. Go is 
similarly defined as a single variable function. The explicit expressions for these form 
factors are presented later in the paper. 

Although the entire effective action is a non-local object, for theories such as QED 
with a massive electron, it should permit a local expansion in inverse powers of m. A 
crucial feature of the above form of the effective action is that it is indeed manifestly local, 
in the sense that the form factors G n have an expansion in positive powers of the D 2 ^. 
This depends on making the particular choice of basis operators above. 

A very general effective action calculation, which encompasses the more specialised 
results we need, has been given some time ago by Barvinsky, Gusev, Zhytnikov and Vilko- 



visky (BGZV) in ref. [T^], building on methods developed earlier by two of the authors 
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n,15,16]. Closely related results have also been obtained by Avramidi | I7| . The present 



paper is essentially a translation and specialisation of the BGZV result to the O(RFF) 
terms in the QED effective action. However, the BGZV result is presented in an appar- 
ently non-local form and the translation into a manifestly local effective action through an 
appropriate choice of basis operators involves a number of subtle manoeuvres. Together 
with the complexity of the manipulations on the form factors, we feel this justifies our in- 
dependent presentation, in which we explain the delicate technical points in the translation 
and quote explicit results for the form factors in the local effective action (|1.11|) . 

In the next section, we introduce the BGZV effective action and show how it is recast 
in local form and reduced to the O(RFF) action we are looking for. Then in section 3, 
we reconsider the 0(F 2 ) terms and rewrite them in a convenient form consistent with the 
Drummond-Hathrell action. Section 4 contains a summary of our final result, with explicit 
algebraic expressions for the form factors collected in an appendix. A numerical analysis 
of the form factors is given in section 5. The paper is presented so that the final results 
for the local effective action and the form factors are entirely self-contained in section 4 
and the appendix. 



2. The BGZV Effective Action 



Barvinsky et al. evaluate the heat kernel Q(x,x';s) (denoted by K in refs.|13|,|T5| 
corresponding to the generic second-order elliptic differential operator 



H = ni + P - \R% U = g^D^D v (2.1) 

D 



acting on small fluctuations 8(j) of an arbitrary set of fields (p (x). The hat notation 
denotes matrices acting on the vector space of the 5(j) A , i.e. P = P A b, 1 = 5 A b, 
etc. and the matrix trace is denoted tr. The Euclidean metric g^ v describes the background 
spacetime, with gravitational Ricci and Riemann curvatures R^ v and R^ v \ p respectively. 
A further 'curvature' is defined by the commutators of the covariant derivatives acting on 
the fields 5<p A : 

(D^D, - D v D^6cf> A = n A B ^54> B (2-2) 

and we denote 71^ = 1Z A b^w The 'potential' P is an arbitrary matrix. BGZV introduce 
the collective notation 9\ for the full set of generalised curvatures: 

« = {R^x P , P] (2.3) 

and calculate the heat kernel Q up to Offi 3 ). 
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P = iea^F^ - —Rl (2.4) 



The first step in extracting the QED effective action from the BGZV results is to 
identify the curvatures 9^ by comparing the differential operator V of eqs. (|1.8| ),( |1.10| ) with 
the standard form H. In this case, the field <p A is the electron spinor field ip a (x), so we 
identify A, B as Dirac spinor indices a, (5 and identify 

1 

12" 

The (gauge and gravitational) covariant derivative acting on spinors is 

D„ = + ^a ab uj abtl + ieA^ (2.5) 
where w a 6 M is the spin connection 

^V = e^(d M e%-r^e a A ) (2.6) 

with e a M the vierbein and a ab = |[7 a ,7 & ]. Evaluating the commutator, we find 

D„] = ieF^l + -a ab R abfl „ (2.7) 



and so, comparing with eq.( |2.2|) , we identify 

%v = ieF^l + X -a xp R^ p (2.8) 

Having established this identification between the BGZV generalised curvatures and 
the QED field strength and gravitational curvature tensors, we can now summarise their 
result for Tr Q(x,x';s). Here, Tr denotes the full functional trace, i.e. Tr Q(x,x';s) = 
J dxy/g tr G(x, x r ;s)\ ,_ . The result is [[13 1 : 



11 _ 

TrQ = (4^p^ / dx ^ tr 



1 + sP 
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i=l 
11 

25 

+ s 4 ^(-50(1), -sD (2) , -SD( 3 )) ^1^2^ 3 (^) 
i=12 
28 

i=26 

+ s 6 F 29 (- S n (1) ,- S D (2) ,- S n (3) ) ^^2^3(29) + 0(* 4 



(2.9) 



Here, IHiD^W denote independent basis operators formed from terms of second or- 
der in the generalised curvatures ( [2.3| ). Similarly for the independent third order terms 
DlitRffis (*)• Notice that these appear with different powers of s, determined by their 
dimension. The form factors and Fi are functions of the operators etc. with 

the suffix indicating which curvature term is acted upon. Thus, for example, since 
9^1^2(5) = ii-iij,v'li-2^ v ' 1 the corresponding contribution to Tr Q is 

J dxy/g n^M-sO)^ (2.10) 

We will return to the detailed expressions for the form factors later. First, we consider 
the basis of 5 0(9\ 2 ) and 29 0(9l 3 ) operators and, with the identifications in eqs. (|2.4|) and 
( |2.8|) , pick out those from the full list which produce the terms of 0(F 2 ) and O(RFF) 
relevant for our QED effective action. The remaining terms involving purely gravitational 
curvatures of 0(R 2 ) and 0(R 3 ) are neglected. In BGZV notation, the relevant basis 
operators are as follows: 

tr 9^2(4) = trAA = \e 2 F pu F pv tri + 0(R 2 ) 

tr 9ti9t 2 (5) = tr Hx^H^ = -e 2 F^F^ tri + 0{R 2 ) 
tr9*i ^2^3(1) = tr P 1 P 2 P 3 = -^(RF^F^ + F pu RF» v + F pu F pv R) tri + 0(R 3 ) 

triRi ^2^3(3) = tr£ lM ,£fP 3 = (^F^F^R+^R^xpF^F^ 

+ \e 2 F pv R pvXp F Xp } tri + 0(R 3 ) 

tr JRi«2JH 3 (6) = trAAi?3 = \e 2 F^ v F pu R tri + 0(P 3 ) 

tr 9^2^3(7) = tr R{k 2vLV 1l pv = -e 2 RF^ v F pv tri + 0(i? 3 ) 

tr ^^2^3(8) = tr Ri^TZ^iZ^x = -e 2 R^F pX F v x tri + 0(R 3 ) 

tv^ 2 fR 3 (U) = tr D^ X D„K V 2X P 3 = {^e 2 D^D v F v xR+ X -e 2 D^R pa ^D v F v x FP a 

+ ^e 2 D^ x D u R pa » x F^ tri + 0(R 3 ) 

X -e 2 R pv D p D v F Xp F Xp tri + 0(R 3 ) 

= -e 2 R PLU D x F xp D p F pu tri + 0{R 3 ) 
= -e 2 R pv D p F xp D u F Xp tri + 0{R 3 ) 
= -e 2 RD VL F pX D v F v x tri + 0(R 3 ) 
= -e 2 R pv D p D x F Xp F pv tri + 0(R 3 ) 

(2.11) 



tr 91^2^3(17) 


= tr 


R 1[1V D»D V P 2 P 3 = 


tr 9*i9t 2 9t 3 (18) 


= tr 


R lfliy D x n Xp D p TZ pu 


tr 9^i9^ 2 9^3(19) 


= tr 


R lpu D p 'R Xp D»'R 3Xp 


tr 9^i9l 2 9^3(20) 


= tr 


R l D^n px D u n 3 ^ x 


tr 9*i9^9t 3 (21) 


= tr 


R lpu D p D x n 2Xp n pi/ 
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Notice that for clarity we have omitted the subscripts 1,2,3 on the terms in the right-hand 
column, but must remember that the order is still significant since they are to be acted on 
by the form factors. 1 is the unit Dirac matrix, so in four dimensions trl = 4. 



Collecting these terms and substituting into eq . (pT9|) then gives the complete 0(F 2 ) 
and O(RFF) terms in the BGZV effective action. However, as it stands this is expressed 
in non-local form, whereas we require a manifestly local action. As explained in ref. Jl5[] , 
the key is to reintroduce the Riemann tensor into the basis of operators. (Notice that none 
of the above set of terms in the form tr D^i^^ involves the uncontracted Riemann tensor 
Rnv\pi which only enters in the right-hand column when we expand TZ^ and P according 
to eqs.(|2T4T),(|2.8|). This is a feature of the covariant perturbation theory construction 
developed by Barvinsky and Vilkovisky.) We therefore introduce the further operator 



tr Ri^ uXp 1Z p ] u iZ Xp 



-e 2 R^ uXp F piJ F Xp tri + 0(R 3 



(2.12) 



The Riemann tensor admits the following non-local re-expression in terms of the Ricci 
tensor, up to 0(R 2 ): 

R^\ P = j=j {d^xRup — P>vP>\Rp P — DpDpRvX + D^DpRfixj + 0(R 2 ) (2.13) 



This is proved [15 1 by differentiating the Bianchi identity 



R 



pu[\p;a] 



(2.14) 



to obtain an equation for DR pu \ pi which is solved iteratively in powers of curvature using 
the Green function 1/D to obtain (|2.13| ). 

We can then rewrite the four BGZV basis operators which turn out to have non-local 
form factors as: 



/ 



dx^/g tr 



sFsR^U^x + s 2 F 18 R^D x n x ^D p np» 

+ s 2 F 19 R^D^TZ Xp D"TZ Xp + s 2 F 21 R^D p D x n Xp TZ p " 



dxyfg tr 



t pv 



sF 8 R^ti» x K u x + s 2 F 18 R^D x n x »D p W 
+ s 2 F 19 R^D^TZ Xp D"TZ Xp + s 2 F 21 R^D p D x TZ Xp n p " 

+ sF R^ Xp iz^n Xp 



(2.15) 
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with the following relation between the form factors: 



F 8 (xi, x 2 , x 3 ) 
F 18 (x 1 ,x 2 ,x 3 ) 
F 19 (x 1 ,x 2 ,x 3 ) 

F 21 (x 1 ,X 2 ,X 3 ) 



X 2 + X 3 
X\ 



F 8 (xi,x 2 ,x 3 ) + 2(l + 

4 „ 

F 18 (xuX2,x 3 ) - — F (x 1 ,x 2 ,x 3 ,) 

X\ 

2 ~ 

F 1 q(x 1 ,x 2 ,x 3 ) + —F (x 1 ,x 2 ,x 3 ) 

8 

F 21 (x 1 ,x 2l x 3 ) - — F (x 1 ,x 2 ,x 3 ) 

X\ 



F (xi,x 2 ,x 3 ) 



(2.16) 



where we have used the notation Xi = — Ouy (Notice that BGZV use £j = — sD(i) in 
discussing the form factors.) 

The new form factor Fq(x\, x 2 , x 3 ) is defined so as to remove the singular terms in 
^8, 18, 19, 21 leaving the Fg, 18,19,21 form factors regular as Xi,x 2 ,x 3 — > 0. Of course this 
under-determines Fq, but any ambiguity only amounts to a reshuffling of contributions 
amongst the various terms in the final expression for the effective action. 

In the local basis, we therefore find that all the 0(F 2 ) and O(RFF) terms in the heat 
kernel are contained in the following BGZV terms (again suppressing the arguments of the 
form factors, i.e. /j = /j(-sD( 2 )) and = F^-sD^, -sD( 2 ), -sD( 3 ))): 

dx^tr /4A-P2 + hnx^n^ 

+ sF 1 P 1 P 2 P 3 + sF 3 £ lM „£fP 3 + sF 6 p!P 2 R 3 + sF 7 R{k 2 ^Tif 

+ s 2 F 14 D^kfD v 1l v 2X P 3 + s 2 F 17 R lpv D»D v P 2 P 3 + s 2 F 20 RxD^VJ^DJR 

+ sF 8 Ri^TZ^ 7Z 3 \ + s 2 Fi8 R\fj, v D\R, 2 lJ l D p Ti p 3 

+ s 2 F 19 R 1 ^D^ p D'^R 3Xp + s 2 F 21 R 1 ^D x n 2Xp n pu 

+ sF R^xpH^K^ 

(2.17) 

We can now re-express this explicitly in terms of the 0(F 2 ) and O(RFF) basis op- 
erators of eq. (|l.ll|) using the set of relations (|2.11| ). We also need to use the following 
identity, which is proved by integrating by parts and using the Bianchi identity for F pv : 

J dx^ D x R VLvXp D a F ap F^ v 

= - J dx^ ( R pvXp D (J F° p D x F^ + l -R^ Xp UF^F Xp + 0(R 2 ) ) 

(2.18) 



V 

3 A 
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This gives: 



Tr g 



(An) 



rtrl 



dx^/g 



F» v h 

+ s fa RF^F"" + h 2 R^F^F'x + h 3 R^ Xp F^F x ^ 

+ s 2 (h 4 RD p F» x D v F v x 

+ h 5 R pv D x F x »D p F^ + h 6 R pv D»F x »D v F Xp 
+ h 7 R pv D»D v F Xp F Xp + h 8 R pv D p D x F Xp F pv 

+ h 9 R puXp D a F^D x F^ ) 

(2.19) 



where ho = ho(— sD) and hi = hi(—s\3^, — sD(2), — sD(3)), i > 1. The form factors are 
related to the BGZV definitions as follows (where we assume appropriate symmetrisations 
on x 1 ,x 2 ,x 3 ): 



h {x) = -^h(x) + f 5 (x) 
h 1 (x 1 ,x 2 ,x 3 ) = -F 1 (x 1 ,x 2 ,x 3 ) - —F 3 (x 2 ,x 3 ,x 1 ) 

- -^F 6 (x 2 ,x 3 ,x 1 ) + F 7 (x 1 ,x 2 ,x 3 ) 

h 2 (x!,x 2 ,x 3 ) = F 8 (x 1 ,x 2 ,x 3 ) 

1 1 

h 3 (x!,x 2 ,x 3 ) = --F 3 (x-l,x 2 ,x 3 ) - -F 14 (x 1 ,x 2 ,x 3 )(x 2 + x 3 ) 

+ F (x 1 ,x 2 ,x 3 ) 

h 4 (x 1 ,x 2 ,x 3 ) = -—F 14 (x 2 ,x 3 ,x 1 ) + F 2 o(x 1 ,X2,x 3l ) 

h 5 (x!,x 2 ,x 3 ) = F 18 (x 1 ,x 2 ,x 3 ) 
he(x!,x 2 ,x 3 ) = F ig (x 1 ,x 2 ,x 3 ) 

h 7 (x!,x 2 ,x 3 ) = --F 17 (x 1 ,x 2 ,x 3 ) 

h 8 (x!,x 2 ,x 3 ) = F 21 (x 1 ,x 2 ,x 3 ) 

hQ(x!,x 2 ,x 3 ) = ^F 14 (xx,x 2 ,x 3 ) 

(2.20) 
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3. Second Order Terms 



While eg. ( |2. is a perfectly satisfactory, local expression for Tr Q, it is not the most 
convenient form for many applications, e.g. photon propagation. Nor is it the form which 
permits the simplest comparison with standard Schwinger-deWitt results. The key step in 
transforming to this standard form is to rewrite the second-order contributions of the type 
/ dxyjg F^UF^ as / dx^/g D^F^ X D„F V A , plus higher order terms of O(RFF) which 
of course modify the corresponding form factors. This second choice for the 0(F 2 ) basis 
operator also matches the term in the Drummond-Hathrell action (|1.2|) and turns out to 
be the more useful representation. 

This transformation can be made already at the level of the BGZV curvatures using 
the identity |T5| , |T3| : 



tr J dx^ % v HRJ lv = tr J dx^ ( - 2D^ x D v n v x + AU^K X P K P ^ 

+ 2R llv n lJ ' x n v x - R^ Xp n^n Xp 



(3.1) 

In terms of the electromagnetic field strength and gravitational curvature, this is equivalent 
to 



J dx^ F^UF^ = J dx^ ( - 2D P F» X D V F V X 

+ 2R^ X F V X - R^ Xp F^F Xp + 0(R 2 



(3.2) 



The proof is straightforward, relying on the Bianchi identity for F^ v and the cyclic sym- 
metry of the Riemann tensor. 

We now need to extend this identity to include the full form factors. The first step is 
to separate off the constant term in ho(— □) and define 

j dx,/g F^h (-O)F^ = MO) / dxy/g F^F^ + J dx^g i^-DpF^ (3.3) 

The first term on the r.h.s., where h(0) = —1/6, eventually gives rise to a divergent 
contribution to the effective action which is removed by the conventional electromagnetic 
field renormalisation. The remaining term simplifies in two stages. First, write 

J dxy/g F^h(-n)BF^ = J dx^ ( - 2D p F px h{-U)D u F u x + 2R pv F px h{-U)F v x 

-R pvXp F^h{-U)F Xp + 2D fX F^ x [D l/ ,h(-D)]F\^j 

(3.4) 
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The second and third terms are already in the required form of O(RFF) basis opera- 
tors. Neglecting higher-order terms in the curvature simplifies the commutator, and with 
repeated use of the Bianchi identity for the Riemann tensor we eventually find 

J dx^ F^h{-U)UF^ = j dx^ (-2D p F pX h(-D)D„F\ 

+ 2R^(h(-n) + nti(-n))F» x F»\ - R puXp h{-u)F^F Xp 

-2R pu h\-U)D x F Xp D p F pv + AR pu h\-U)D p D x F Xp F pv 
+ 2R puXp h\-U)D (J F ap D x F pv + 0(R 2 FF)^j 

(3.5) 

The transformation between the alternative choices for the 0(F 2 ) basis operator therefore 
affects the O(RFF) form factors h2,h 3 ,h B , hg, and hg , in the notation of ( |2.19| ). 



Tr g 



With this transformation, we can rewrite Tr Q into the preferred form: 

-\f uv F^ + s D^F px g D V F V , 



(4tt) : 



rtrl 



dxy/g 



6 



A 



+ 8 (#i RF^F"" + g 2 R p „F pX F v x + g 3 R puXp F pv F Xp ^ 

+ s 2 (#4 RD p F pX D v F v x 

+ g 5 R pv D x F xp D p F pv + g G R pv D p F Xp D v F Xp 



+ g 7 R pv D p D v F Xp F Xp + g 8 R pu D p D x F Xp F 



+ g 9 R pvXp D a F° p D x F pv 



pv 



(3.6) 

with go = go(— sDi) and gi = gi(— sDm, — sD( 2 ), — sd(3)), % > 1 as before. The new form 
factors are related to the initial ones by: 

2/1 1 

go{x) = -2h(x) = - {--U(x) + fe(x) + - 

x V 2 b 

91,4,6,7 = ^1,4,6,7 

g2 = h 2 + h(x 2 ) + x 2 h'(x 2 ) + h(x 3 ) + x 3 h'(x 3 ) 

93 = h 3 - -(h(x 2 ) + h(x 3 )) 

9b = h B - (h'(x 2 ) + h'(x 3 )) 

9s = h 8 -4ti(x 2 ) 

g 9 = h 9 + 2h'(x 2 ) 



(3.7) 



where g^ = gi(xi, x 2 , x 3 ), i > 1 and we have symmetrised where appropriate. 



11 



4. The Local Effective Action 



Collecting these results and evaluating the effective action from 

1 f°° r/s 2 

T = T - - -e~ m s Trg(x,x>-s) 

2 Jo s 

we find the renormalised expression: 



(4.1) 



r = J dx^/g 



+ G[ RF^F^ + (T 2 R^ X F V X + R^xpF^F^ 

+ -^(gIrd p f» x d u f» x 

+ Gl R pv D x F x ^D p F^ + Gl R^D»F x ?D v F Xp 
+ G? R^D v F x PF Xp + Gl R pv D»D x F Xp F^ 

+ Gl R pvXp D a F^D x F^ ) 



(4.2) 



where the form factors GC = G n ( — 

of the g n , viz. 



2 5 



~t- 5 — ^t") are defined as proper time integrals 



G n (x 1 ,x 2 ,x 3 ) = -\— ( —e s s p g n (sx 1 ,sx 2 ,sx 3 ) 

2 7T Jo s 



(4.3) 



where p = 1 for n = 0, ... 3 and p = 2 for n = 4, ... 9. Note that we have rescaled s 
here so that it is dimensionless, while now Xi = -□ (i) /m 2 . The divergent contribution 
to the F^F^ term has been evaluated using dimensional regularisation and cancelled by 

the standard electromagnetic field renormalisation. The residual finite 0(a) contribution 

i 2 

gives Z = 1 + 77— In ^ for a suitable choice of RG scale a. 

The form factors are given in terms of the BGZV definitions by the following set of 
relations: 



9o{x) 



-2h(x) = \{-\Sa{x) + h{x) + 1) 



gi(x 1 ,x 2 ,x 3 ) = ^F 1 (x l ,x 2 ,x 3 ) - -^F 3 (x 2 ,x 3 ,x 1 ) 

-2^6(^2,^3,^1) + F 7 (x 1 ,x 2 ,x 3 ) 
g 2 (x 1 ,x 2 ,x 3 ) = F 8 (x 1 ,x 2 ,x 3 ) + h(x 2 ) + x 2 h'(x 2 ) + h(x 3 ) + x 3 h'(x 3 ) 

12 



g 3 {x x ,x 2 ,x 3 ) = -~F 3 (xi,x 2 iX 3 ) - ~Fi 4 (xi,x 2 ,x 3 )(x 2 +x 3 ) 

+ F (xi,x 2 ,x 3 ) - -(h(x 2 ) + h(x 3 )) 

5r 4 (xi,x 2 ,x 3 ) = -^-Fi 4 (x 2 ,x 3 ,xi) + F 2 o(xi,X2,X 3 ) 

5r 5 (xi,x 2 ,x 3 ) = Fi 8 (xi,x 2 ,x 3 ) - (h'(x 2 ) + h'(x 3 )) 

fir 6 (xi, X 2 , X 3 ) = Fiq(xi,X2,X 3 ) 

fir 7 (xi,X 2 ,X 3 ) = -^Fi 7 (xi,X 2 ,X 3 ) 

g8(x 1 ,x 2 ,x 3 ) = F 2 i(xi, x 2 , x 3 ) - Ah'{x 2 ) 

fir 9 (xi,x 2 ,x 3 ) = Jf m (xi,x 2 ,x 3 ) + 2h'(x 2 ) 



(4.4) 



where h'(x) = 4^ and the are given in terms of the BGZV Fi by eq. (|2.16| ). Again we 
assume that the BGZV form factors have been appropriately symmetrised to reflect the 
symmetries of the basis operators on which they act. 

The BGZV form factors themselves are known algebraic expressions which involve the 
following integrals []X5|JT^] : 

f(x) = [ da e"^ 1 -^ (4.5) 

and 



o 



F(xi,x 2 ,x 3 ) = f rf 3 a 5(1 -«i -a 2 -a 3 ) e- (QlQ2X3+Q2a3a;i+Q3aiX2) (4.6) 

Jot>0 

The second-order form factors are easily simplified. Using the BGZV expressions, 

h{x) = \f{x) hix) = ~(f{x)-l) (4.7) 

we have 

Using the small x expansion of f(x), viz. 

f( x ) = 1 - -x+ — x 2 - — x 3 + . . . (4.9) 

Jy 1 6 60 840 K ' 

we find 

g (x) = x + . . . (4.10) 

yuv ) 15 14Q v ) 
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The absence of inverse powers of x confirms that this is indeed a local form factor. 

The expansions for the third-order form factors F{(xi, x 2 , £3) are extremely compli- 
cated. They can be read off from the results given in ref. []13[ and we will not quote them in 
full here. Instead, we focus on a special case of particular interest. In most applications, 
it is reasonable to assume that the derivatives of the gravitational curvatures are governed 
by the same scale as the curvatures themselves, i.e. if = O(j^) then ^§ = O (-£§■). 
In that case, terms with derivatives of curvatures are suppressed by the same parameter 
governing the expansion in increasing powers of curvature, and should be neglected. We 
would then only be interested in the form factors with the first argument set to zero, i.e. 

fir n (0, £ 2 ,£ 3 )- 

We have calculated the corresponding form factors Fi(0, £2, £3) an d Fi(0, x 2 , £3) ap- 
pearing in eq. (|4.4|) , identifying the new form factor Fq(0, x 2 , £3) required for the transfor- 
mation to the local basis. The resulting lengthy expressions are collected in the appendix. 

These results must pass three crucial consistency checks. First, the transformation 
from the non-local basis to the local one via the introduction of Fq(xi, x 2 , £3) as described 
in eq.( |2.16| ) must work, i.e. the l/xi singularities in each of the form factors F 8j 18,19,21 
must be removed by the same function Fq. Second, all the resulting form factors must be 
entirely local, in the sense that they admit Taylor expansions with no inverse powers of 
£1,2, 3- Finally, the final form factors <7o(0) and gr n (0, 0, 0), n = 1,2,3 must reproduce the 
Schwinger-de Witt coefficients in the Drummond-Hathrell effective action. In algebraic 
terms, each of these tests is highly non-trivial, but as can be verified from the explicit 
results in the appendix, each is indeed satisfied. 



The £1,2,3 - ► limit of the form factors in the effective action can be found from the 
results quoted in the appendix, using G n (0,0,0) = -if £ n (0,0,0). From the definitions 
( [PP we find: 

Go(0) = ~- Gi(0,0,0) = G 2 (0,0,0) = ~ 

30 TV 144 TV 360 TV 

G 3 (0,0,0) = -~ G4(0,0,0) = -L- G 5 (0,0,0) = -L- 

360 tv 3360 7r 315 tv 

G 6 (0,0,0) = -—- G 7 (0,0,0) = — - G 8 (0,0,0) ' 



G 9 (0,0,0) 



2520 7T v ' ' ' 720 7T v 315 tt 

5 a 



1008 tv 



(4.11) 



and we confirm that Go(0), (2^2,3(0,0,0) correctly reproduce the coefficients d, a, 6, c re- 
spectively in the Drummond-Hathrell action ( |1 . 2|) . 
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5. Numerical Analysis of the Form Factors 



In order to get a better feeling for the behaviour of the form factors, we have evaluated 
them numerically. We present some of the results in graphical form in this section. The 
general features of the different form factors are quite similar, so here we just give the 
analysis for the form factors relevant for Ricci flat spaces, i.e. #3(0, £2, £3) and #9(0, X2, £3). 

Collecting formulae from the appendix and substituting into eq.(fO|) we find 



fir 3 (0, X 2 ,X 3 ) 



F(0, x 2 , x 3)-^- 4 + (l + -(x 2 + x 3 ))(x 2 + x 3 ) 



+ /(* 3 )i 



3^2 - x 3 1 



3x 



3 ^2 1 / V , 

— + Z (X3-X 2 )(1- 



+ 



A L 



1 /X3 X2 

-1+77 1 h x 2 +x 3 

2 Vx 2 x 3 



2(x 2 + x 3 ) 2 
A 

2(x 2 + x 3 ) 2 
A 
1 



1/1 1 

12Vx 2 x 3 



1/1 1 

4\Xn xi 



(5.1) 



and 



5r 9 (0,X 2 ,X 3 ) 



F(0, x 2 , x 3 ) — [4 + x 2 + x 3 ] 



2x| 



2 1 



Xr 



r 2 1 /ii 

+ /( x 3)[-^2(^ -X^)J + ^( X2 K^2 + ^T 



2 



1 2 



(5.2) 



The behaviour of these functions is illustrated in the following plots of #3(0,X2,X3) 
and (79(0, X2, X3). The values at the origin are ^3 (0,0,0) = j|q and ^9 (0,0,0) = -^g, as 
already quoted in eq . (f4 . 1 1| ) . Notice that along the diagonals X2 = x 3 , both functions tend 
asymptotically to zero. This is consistent with the vanishing of the BGZV form factors 
as □ — > 00 as can be seen from the definitions in eqs.( [4.5|) and Q4.6|) . However, if one 
argument is set to zero, then the functions may tend to a finite limit. It is easy to check 
that this behaviour is entirely determined by the h(x) and h'(x) factors in eq.( |4.4j ). These 
limits are #3(0, 00, 0) = #3(0, 0, 00) = #9(0, 00, 0) = and #9(0, 0, 00) = ^ as can be 
seen in the figures below. 
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10 



100 



Fig. 2 3D plots of 93(0, X2, £3) over different ranges. 

. 014r 

0.012- _______ ■ ' 

0.01- ^ 
. 008 - 
0.006 / 
. 004 - 
. 002 

20 40 60 80 100 20 40 60 80 100 

Fig. 3 Graphs of 93(0, X2, £3) along the X2 or X3 axes (left) and the diagonal X2 = 2:3 (right). 




Fig. 4 3D plots of 99(0, X2, £3) over different ranges. 



.007 
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.005 
.004 
.003 
.002 
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Fig. 5 Graphs of 39(0, X2, X3) along the X2 axis (top left), the diagonal X2 = £3 (top right) and the X3 
axis (lower). 



This completes our numerical discussion of the form factors. We confirm that they 
are regular at Xi — as required for their role in the local effective action. Their values at 
X{ = match the well-known results for the Drummond-Hathrell effective action. They are 
smooth functions with well- under stood asymptotic behaviour. The local effective action 
( |4.2| ) for photon-gravity interactions induced by vacuum polarisation is therefore under 
complete algebraic and numerical control to all orders in the derivative expansion. The 
precise behaviour of the form factors may now be used to address a variety of physical ap- 
plications, notably the important issue of dispersion for photon propagation in background 
gravitational fields (1^ . 
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Appendix A. 



In this appendix, we quote explicit formulae for the third order form factors appearing 
in the heat kernel in the limit discussed in section 4. 

These are expressed in terms of the two basic integrals 

f{x) = [ da e - ail ~ a)x (A.l) 
Jo 

and ^ 

F(0,x 2 ,x 3 ) = f da - f e - Q(1 - Q)x3 _ e -«(i-«)*a > ) (A.2) 

x 2 - x 3 J a \ J 

We will at times use the further limit 



F(0,x,0) = - I da-(l-e- a{1 - a)x 
x Jo ol \ 



11 1 2 1 , 

x H x x + 

2 24 360 6720 



The final expressions also involve 



X + . 



1 


1 


6x 


2x 2 


1 


1 






30 + 


280 



(A.3) 



d 

F'(0,x 2 ,x 3 ) = —F(x 1 ,x 2 ,x 3 ) (A.4) 

OX\ xi=0 

which has a small x expansion 

F'(0,x,0) = -— + — x —x 2 + ... (A.5) 

v ' ' ' 24 360 6720 V ; 

We also define A = (x 2 — x 3 ) 2 . 

It is convenient to recall here the definition of h(x) from eqs.(3]7),(0), viz. 



(A.6) 

Beginning with the form factors which are not changed in the transformation to the 
local basis, we need the following contributions to gi(0, x 2 , x 3 ). Note that these have been 
appropriately symmetrised. 

F 1 (0,x 2 ,x 3 ) = ^F(0,x 2 ,x 3 ) (A.7) 
18 



2 2 

F 3 (x 2 ,x 3 ,0) = F(0, x 2 ,x 3 ) -^-x 2 x 3 - — (x 2 +£3) 



/ \ 4 1 

\J{x 2 ) ~ f{x 3 )j-^X 2 X 3 {x 2 - X 3 ) + — [X2 + X3) 



(A.8) 



1 1 

F 6 (x 2 ,x 3 ,0) = F(0,x 2 ,x 3 ) - — (x 2 + Ax 2 x 3 + x\) - — (x 2 + x 3 ) 



f(x 2 )-f(x 3 ) 



-^(x 2 -x 3 )(xj + 6x 2 x 3 +xl) + ^(x 2 + x 3 ) 

(A-9) 



F r (0,x 2 ,x 3 ) = F(0,x 2 ,x 3 ) 



^^x 2 x 3 {x 2 2 + 4x 2 x 3 + x\) + 7^(>2 + x 3 )(x\ + 34x 2 x 3 + x\) 



+ (f(0,x 2 ,x 3 ) - ^)-^( x 2 + 10x 2 x 3 + x\) - ~^{x 2 + x 3 )x 2 x 3 



+ 



+ 



1 / / x 2 1 \ J_ X2 /a _|_ 44 a .3 a , 3 _ io£ 2 a; 3 - 36x 2 a; 3 + xi) + x 2 <-> x 3 

2 V jo / A* 5 



^ ( f ^x 3 " ) a^ 3 ^ 2 " 36x3a; 2 - !0^^2 + 44x^2 + x£) + x 2 



^3 



(A.10) 



For (74(0, x 2 , £3) and <7g(0, x 2 , x 3 ), we require the following: 



Fi 4 (x 2 , x 3 ,0) = F 14 (0,x 2 ,x 3 ) 

2 / 1 \ 8 

= F(0,x 2 ,x 3 )-(x 2 +x 3 ) + (F(0,x 2 ,x 3 ) - -J- 

+ (/(x 2 )-/(x 3 ))^(^-^) 



(A.ll) 



F 20 (0,a; 2 ,x 3 ) = F(0,a; 2 ,a; 3 ) 



1 2 

(x 2 + x 3 )(x2 + 4x 2 x 3 + £ 3 ) - ——^(17x1 + 38x 2 x 3 + 17x 3 ) 



3A 2V /w 6J 3A 2 

2 

/Mg^g (>2 - X 3 ){X 2 + X 3 )(x 2 + Ax 2 X 3 + X 3 ) + X 2 <-> X 3 

ffix?) 1\ _L( X2 _ X3 )(21x 2 + 45x|x 3 + 15x 2 xi - x 3 ) + x 2 <-> x 3 

V x 2 / A d 

(A.12) 
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For #3(0, X21 x 3 ) we need 

F 3 (0,x 2 ,x 3 ) = -(f(x 2 ) - f(x 3 )j — (x 2 - x 3 ) 
while #7(0, X2, x 3 ) involves 

F 17 (0,X 2 ,X 3 ) = -F(0,X 2 ,X 3 )^ + [f(x 2 ) + f(x 3 ] 



1 

A 



(A.13) 



(A.14) 



We have explicitly checked that these form factors are indeed local by evaluating the 
limits Fi(0, x, 0) and verifying the cancellation of all the terms with inverse powers of x. 

The remaining form factors involve the transformation ( |2.16| ) between the original 
non-local BGZV form factors and the new local ones. The first step is to identify the new 
form factor Fq(0, x 2: x 3 ). We find this from any of the identities in eqs.( |2.16| ). For example 
(remembering that we understand Fig, to be symmetrised on x 2 ,x 3 and using the BGZV 
expression |L3|]), 

F (0,x 2 ,x 3 ) = -- lim Xi F 18 (xx,x 2 ,x 3 ) (A. 15) 

All the identities produce the same result for Fo(0, x 2 , X3), verifying the consistency of our 
method. We find 

-)- 

2J A 

— (3x 2 -x 3 ) - (f(x 3 )-l)— — 
x 2 \ J 2A x 3 

(A.16) 

The remaining form factors are then: 



F (0, x 2 ,x 3 ) = -\F(0,x 2 ,x 3 ) 

+ ( /(X2) - 1 )2ki (3X2 ~ X3) ~ ( /(X3) - 1 )2Z^ (X2 - 3XJ 



F 18 (0, x 2 ,x 3 ) 



20 

F(0,x 2 ,x 3 )^2(x 2 + x 3 



2 

1\ 80 
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2 j A2 \X2 + x 3 ) + — F'(0,x 2 ,x 3 ) 



Fi 9 (0,x 2 ,x 3 ) = 



+ (f(0,x 2 ,x 3 ) 

+ 2 (^-^ -^-^(x 2 -x 3 )(l7xl + 45xjx 3 + 19x 2 x 2 3 -x 3 3 ) + x 2 <-> x 3 

(A.17) " 

- F(0,x 2 ,x 3 )^(x2 + 3x 2 x 3 + x\) 



+ 



(F(0, x 2 ,x 3 ) - ^)^( x 2 +x 3 ) - —F'(0,x 2 ,x 3 ) 

f(x2) ~ 



X 2 



--^xj(x 2 - x 3 )(3x 2 + 7x 3 ) 
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+ -^jx 2 (x 2 + x 3 ) (x\ - X2X3 + 6x3X3 - 4x 2 x 2 + X3) ] + x 2 



x 3 



(A.18) 
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1 6 

i(0,x 2 ,x 3 ) = F(0,x 2 ,x 3 )-^(x 2 + x 3 )(x 2 + 4x 3 ) 

1\ 64 



+ (f(0,x 2 ,x 3 ) - i)^(aJ2 + 4x 3 ) + ^F'(0, x 2 ,x 3 



and finally 

F 8 (0, x 2 ,x 3 ) 



+ I ~)(^3"( X 2-^3) (9^2 + 21^3 + ll^ 2 xl -dj) - ^X 2 (X 2 +X 3 ; 

f I A^a) — 1^^^^ _ xz){x% - 6x2x3 - 5x§) - ^-x 3 (x 2 + x 3 ; 



X3 



A 2 ' 



(A.19) 



40 

F(0,x 2 ,x 3 )-^x 2 x 3 (x 2 + x 3 ) 

x 1\160 
F(0, x 2 , x 3 ) - - ) ^-x 2 x 3 



+ 



f(x 2 ) - 1 



x 2 



-(x 2 +x 3 )F'(0 
-^3 (^2 _ x 3 )(l7x 2 + 86x 2 x 3 + 72x 2 x 3 - 22x 2 x 3 + 7x 3 ) 



+ ^i( X 2 + ^3) (3X2 + 2X 2 X 3 - X3) J +x 2 ^x 3 



(A.20) 



Again, using the results above, we have explicitly checked that the form factors 
^0,8, is, 19,21 all have regular small x limits, as required for the local operator basis. The 
required numerical cancellations are complicated and provide a very stringent test of the 
algebraic forms quoted. 

Their remaining finite values when all of xi j2j3 — ► reproduce the coefficients in 
the Drummond-Hathrell action. To find these, we need the small x expansions of f(x), 
F(0, x, 0) and F'(0, x, 0) given in eqs. Q ,(gg), (gJ ). We find: 



F o (0,0,0) 

and similarly 

Fi(0,0,0) = i 
6 

F 8 (0,0,0) = 



lim 



1 



A 9 (0,0,0) 



180 
1 

1260 



-— (F(0,x,0)--) + 

F 3 (0,0,0) = 
F 14 (0,0,0) : 
F 20 (0,0,0) 



/(x)-lx 3 t 1 



x 



2x + 12x 



— (A.21) 
72 V ; 



1 

12 
1 



180 



F 6 (0,0,0) = 
F 17 (0,0,0) 



1 



7560 



F21 (0,0,0) 



180 
1 

630 



F 7 (0,0,0) = 
F 18 (0,0,0) 



1 



1260 



(A.20) 
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